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Abstract 

In this article, we extensively develop Carleman estimates for the wave equation and give 
some applications. We focus on the case of an observation of the flux on a part of the boundary 
satisfying the Gamma conditions of Lions. We will then consider two applications. The first one 
deals with the exact controllability problem for the wave equation with potential. Following the 
duality method proposed by Fursikov and Imanuvilov in the context of parabolic equations, we 
propose a constructive method to derive controls that weakly depend on the potentials. The 
second application concerns an inverse problem for the waves that consists in recovering an 
unknown time-independent potential from a single measurement of the flux. In that context, 
our approach does not yield any new stability result, but proposes a constructive algorithm to 
rebuild the potential. In both cases, the main idea is to introduce weighted functionals that 
contain the Carleman weights and then to take advantage of the freedom on the Carleman 
parameters to limit the influences of the potentials. 

Keywords: wave equation, Carleman estimates, controllability, inverse problem, reconstruc- 
tion. 

AMS subject classifications: 93B07, 93C20, 35R30. 

1 Introduction 

The goal of this article is to revisit observability properties in the light of Carleman estimates 
for the wave equation with a potential in a bounded domain. We will present applications of 
the appropriate Carleman estimates in two directions: 

• In control theory on the dependence of the exact controls for waves with respect to the 
potentials; 

• In an inverse problem for the wave equation in which the potential is unknown, where we 
will give a reconstruction algorithm for the potential. 
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1.1 Setting 



Let n be a smooth bounded domain of R", n > 1, and T > 0. We consider the wave equation 
d'iz - Az+pz = g, inrix(0,2r), 

2 = 0, ond9.x (0,2r), (1.1) 

z{Q) = Zo, dtz{Q) — zi, in Q. 

Here, z denotes the amplitude of the waves, p is a potential supposed to be in L°°{Q x (0, 2T)), 
g is a source term for instance in L^{Q, x (0, 2r)) and {zq,zi) are the initial data lying in 
H^{Q,)xL'^{Q,). It is by now well-known that, due to hidden regularity results |24) . under these 
assumptions, the normal derivative of z on the boundary belongs to L^{dQ, x (0, 2T)). 
In this article, we focus on the following observability property: 

Given To C dil, can we determine z solution of Hl.l|) from the knowledge of g, p 
and d^z on To x (0,2T)? 

When p = 0, this question has a positive answer if and only if the Geometric Control 
Condition holds [2 E] for ^1,2T and Fq. Roughly speaking, it asserts that all the rays of 
geometric optics in Jl, which here are simply straight lines reflected on the boundary according 
to Descartes Snell's law, should meet the observation region To at a non-diffractive point in a 
time less than 2T. 

However, other methods exist based on multiplier techniques |24l I23| or on Carleman es- 
timates [141 131) . These methods use stronger geometrical assumptions, and in particular the 
following ones, sometimes referred to as the Gamma-condition of Lions or the multiplier con- 
dition. 

Geometric and time conditions: 

3 xo ^ n, such that To D {x e 8^1, {x - xq) ■ v{x) > 0}, (1.2) 
T > supl2;-a;o|. (1.3) 

The advantage of Carleman estimates on the multiplier techniques is that they allow to 
easily handle potentials in L°°(f2 x (0,2r)) - see e.g. |14l 1311 [TT] . In the applications we have 
in mind and that will be developed hereafter, it will be important to understand the dependence 
of the observability inequalities with respect to the potentials. This precisely explains why the 
path we have chosen hereafter uses Carleman estimates. 

In order to state our results precisely, we shall need several notations. To make them easier, 
instead of working on (0, 2T) as in (|l.ip . we do a translation in time in order to consider: 

dtz-^z+pz = g, innx(-r,r), 

z = G, on on x{-T,T), (1.4) 

z(—T) — Zg'^ , dtz{—T) = z^"^ , in f2. 

Let us define, once for the whole paper, the weight functions we shall consider in Carleman 
estimates. 

Weight functions: Assume that To satisfies (|1.2|l for some xo ^ ^- Let /? G (0, 1), and define, 
for {x,t) G SI X (-r,T), 

'il;{x,t) = \x - xo\'^ ~ I3t'^ + Co, and for A >0, (^(i, t) = e^"'"^"'*', (1.5) 

where Co > is chosen such that i/i > 1 in f2 x {—T,T). 

Note that the weight function ip defined that way depends on /3 £ (0, 1) and A > and shall 
rather be denoted ^pp,x, but these dependences are omitted for simplifying notations. 
We also define, for m > 0, the spaces 

L<m{fl) = {q e L°°{^), ||g||i,oo(n) < m}, 

L^^in X (-T,r)) = {pG L°°(f7 X {-T,T)),\\p\\L^(n>,(-T.T)) < m}. 

The main results we shall use are the following ones: 
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Theorem 1.1. Assume the multiplier condition (|1.2p . the time condition (|1.3p . Let /3 G (0, 1) 
be such that 

sup|a:-xo| < I3T. (1.6) 

T/ien for any m > 0, t/iere exisf A > independent of m, so — so(m) > anci a positive 
constant M — M(m) such that for being defined as in (jl.Sp . for all p G L'^^in x (— T, T)) 
anrf /or a// s > Sq: 

s / e^"^ (la^l^lV^l") dxdt + s^ r [ e^'^lz]^ dxdt 

{\dtz{-T)f + \Vz{-T)f)dx + s'' f e''''^^'''^^\z{-T)fdx (1.7) 

Jn 

<M I fe'^'"^\d'}z~Az+pzfdxdt + Msf [ e'^'^ \d^z\'' dadt, 



T Jn 



for all z e L^{-T,T;Hi{n)) satisfying dfz - Az + pz £ L'^ (n x {-T,T)) and d^z G L^(ro x 
(-T,T)). 

Theorem 1.2. Under the assumptions of Theorem M.lX if z furthermore satisfies z[-,Q) — 
in one also has 



j g2s^{0) 

n 



\dtz{0)fdx<M [ [ e^'"^\dtz~^z+pzfdxdt 
J -T Jn 

+ Ms [ [ e^"^ \dyz\^ dodt. (1.8) 

In particular, if z{-,0) = in ft and q G L'^^{ft), then for all z G {0,T; H^){Q,)) satisfying 
dtz - Az + qz G L'^{Q X (0,r)) and d^z G L'^(To x (0,T)) and for all s > so(m), 

^1/2 f e2''^(")|a^^(o)|2rf^ + 5 r f e^"^ {\dtzf + \Vzf) dxdt + s^ f f e^"^\z\' dxdt 
Jn Jo Jn Jo Jn 

<M [ [ e^"^\dU- ^z + qzl^ dxdt + Ms [ [ e^"^ Id^z]"^ dadt. (1.9) 
Jo Jn Jo Jfo 



Note that the condition ^(-,0) = in H of Theorem 11.21 makes sense for z such that z G 
L^({0,T)-Hf,(D.)) a.nd dtZ-Az+pz G L^{nx{0,T)) since then dtZ = Az-pz + {dt z- Az+pz) 
belongs to L^{(0,T)- H-'^{n)). 

Theorems I1.1H1.2I do not claim particular originality and many of their ingredients are 
already available in the literature, see e.g. [2] where very similar estimates are proved and 
[141 1311 117) for more references. However, to our knowledge, this is the first time that these 
global Carleman estimates are written under that form, which is easier to use to achieve our 
goals. Detailed proofs of Theorem ll.lH1.2l are given in Section [21 

1.2 Applications to controllability 

The idea is to take advantage of the Carleman estimate of Theorem 11.11 to obtain controls 
whose dependence with respect to the potential is weak. To be more precise, we focus on the 
following exact controllability problem: 

Given (y,7^,2/f^) G L^{n) x H-\n), find u G L^(ro x i~T,T)) such that the 
solution y of 

( dh~~^y+Py = 0, mQx{~T,T), 

I y^ulro, ondnx{-T,T), (1.10) 

I yi-T) = yo'^, dty{-T) = yi"^, in Vt 

solves 

y{T) = dty{T) = 0, in Q. (1.11) 
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This exact controllability problem is equivalent to the observability of the system (|1.4|) . 
These two properties are dual one from another, as stated by Lions [24 using the so-called 
Hilbert Uniqueness Method (HUM): the HUM computes the control of minimal L^IVq x 
(— T, T))-norm from the minimization of a quadratic functional whose coercivity is equivalent 
to an observability property for the adjoint system H1.4[) that can be deduced from Theorem ll.il 
under the conditions (|1.2|) - H1.3p . 

Actually, we shall not focus on these HUM controls. Nevertheless, the controls that we 
shall consider below are also computed with a duality argument, based on the "observability" 
inequality (|1.7ll directly. Our approach is strongly inspired by the duality strategy employed 
by Fursikov and Imanuvilov |14) . that has mainly been used for parabolic equations so far. The 
idea is to minimize, for s > 0, the functional 

KsAz)^^[ [ e^'^^ldh- ^z + pz\^ dxdt+]- [ [ e'^'^\d,.zf dadt 
2s J -T Jn 2 J_y Jp(j 

+ {{yo'^,yT^), (2:(-r),at2:(-T)))(i2xH-i)x(HixL2), (i-i2) 

on the trajectories z such that z G L'^{-T,T; Ho{0,)), dtZ-Az + pz G L^(n x {-T,T)) and 
d.z G L^{To X (-r,r)). Here, 

r 

{(yo^,yi'^), (2(7^,2r^))(L2xii-i)x(HixL2) = vo'^zi'^ - (yr^, 2,7^)^-1 X Hi. 

with 

{j/r^,2f7'^)^-ixHi = / V(-Ad)"^yr^ • V2(7^, 
Jn 

where Ad is the Laplace operator with Dirichlet boundary conditions. Note that this functional 
Ks,p depends on s - the parameter chosen in the exponential - and on the potential p. 

Then, according to the Carleman inequality in Theorem 11.11 under the conditions (|1.2p - 
(|1.3p . if, for some m > 0, p G L^^{i} x (— T, T)) and s > so(m), we shall easily show that 
Ks,p is strictly convex and coercive. Ks,p therefore has a unique minimizer, denoted by Z[s,p] 
to underline the dependence with respect to the potential p and the parameter s. Simple 
computations prove that if we set 

Y[s,p] = ^e^'"'{d^ - A+p)Z[s,p] and U[s,p] ^ e^"'d,Z[s,p]lro, (1.13) 

we obtain a solution to the exact controllability problem (|1.10p - (|l.lip - see Theorem 13.11 for 
precise statements and proofs. We are then in position to study the dependence of the controls 
with respect to the potential: 

Theorem 1.3. Assume the conditions (|1.2l) - (ll.3p . 

Letm > andp'^jp^ be two potentials in L'^^{Qx{—T,T)). Given initial data (yg'^ jU^"^) G 
L^{Q) X H~^{Q,), there exists a constant M — M{m) > independent of s such that the 
corresponding controlled trajectories {Y[s,p'^],U[s,p"']) and {Y[s,p''],U[s,p'']) satisfy, for all 
s > so(m), 



s 



T 



^ e~^'"^\Y[s,p'']-Y[s,pY dxdt+ [ [ e~''''^\U[s,p''] - U[s,p'']f dxdt 



s 1^ I e 



{Y[s,p''f + Y[s,p''f)dxdt+ [ [ e-'''"^{U[s.,p''f + U[s,p''f)dxdt 

J -T Jvn 



^ ^IIp" -/lli~(nx(-T,T)) < 2mMs-3/2, (1.14) 

where ip = V'/s.a is chosen so that Theorem \1.1\ holds and so(m) is the parameter given by 
Theorem 

In other words, the relative error between the controlled trajectories {Y[s,p],U[s,p]) decays 
as s~^^^ for potentials lying in L'^,„(Q x (— T, T)). 



Theorem 11.31 states that, as s increases, the control obtained by the minimization of Ks,p 
depends less and less of the potential p. This is of course in complete agreement with the 
results obtained using microlocal analysis, in which the potentials play no role - see e.g. [8]. 
However, to our knowledge, the relations between controls computed for different potentials 
have not been studied so far. 

The proof of Theorem 11.31 is purely variational and comes from the variational character- 
ization of the controls. The main issue is to track the powers of s during the proof. Again, 
this strongly relies on the Carleman estimate of Theorem 11.11 One will find all the details in 
Section [3] 

1.3 Applications to inverse problems 

The idea now is to take advantage of the Carleman estimate of Theorem 11.21 to conceive a 
reconstruction algorithm of the potential from the knowledge of the flux of the solution. To be 
more precise, we focus on the following inverse problem: 

Given the source terms h and hg and the initial data (wo,?i)i), considering the 
solution of 

C d'^W - AW + QW = h, inf7x(0,r), 

I W = ha, ondnx(0,T), (1.15) 

[ W{0) = wo, dtW{0) = wi, inn, 

can we determine the unknown potential Q — Q{x), assumed to depend only on 
X £ n, from the additional knowledge of the flux 

fi = d„W, onr„x(0,T) (1.16) 

of the solution? 
Under the regularity assumption 

W € H'-{{0,T);L°°{n)), (1.17) 

the positivity condition 

3a > such that |«;oj > a in Q (1-18) 

and the multiplier conditions p.2p - p.3p . the results in [2] (and in |29) under more regularity 
hypothesis) state the stability of this inverse problem consisting in flnding the potential Q 
from the measurement of the flux (11.161) . To be more precise, it is proved that, given Q", Q'' £ 
L^jj^{n) and denoting by and WfQ*'] the corresponding solutions to (|1.15[l . there exists 

a positive constant M = M{n, T, xq, m) such that 

^IIQ" - Q^h^n) < \\dtd,.wm - dtd,M[Q'']h2io,T;LHro)) < ^WQ" - Q'\\lh»)- (1-19) 

Nevertheless, this stability result (|1.19|l is not given with a constructive argument that 
could allow to explain how to flnd Q from the knowledge of d^Wlro- We are thus interested 
in deriving an algorithm so that it can eventually be implemented into numerical algorithms. 

The algorithm we shall propose is based on a data assimilation problem that we briefly 
present below. 

Let m > and q e -L<„(n). Let € L^iVo x (0,r)) and g £ L'^{n x (0,r)). We introduce 
the functional 

J,,,[M,g](z) = r f e''''^\dfz^Az + qz-gfdxdt + ^ f e'""^\d,.z ~ fi\^ dadt, (1.20) 
Jo Jn ^ Jo Jfo 

on the trajectories z such that z G L^(0, T; J/q (f7)), dtZ - Az + qz e L^{n x (0,r)), d„z G 
L'^iVo X (0, T)) and z(-,0) = in fi. Remark that z(-,0) makes sense in H~^{n) since 
dfZ = Az ~ qz + {dtZ — Az + qz) G L^{Q, T; H~^{0, 1)) under the previous assumptions. 
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Similarly as for ifs.p, one will see that this functional Js,q[iJi.,g\ has a unique minimizer 
Z for s > So (Proposition I4.l| l. We emphasize that Z depends on s and q, but the context 
will make it obvious and we therefore drop these dependences to simplify the notations. More 
importantly, we can study how the minimizer Z depends on g £ L^{Q,x (0,r)) and, similarly 
as in Theorem 11.31 we will prove the following result: 

Theorem 1.4. Assume the multiplier condition (|1.2[l and the time condition (|1.3|) . 

Assume that /i G L^{To x (0,r)) and g", g^ e L^[^ x (0,r)). Let m, > and q G Lf^(r2). 
Let Z-' be the unique minimizer (see Proposition \4.1^ of the junctionals Js.q[^,g-'] for j G {a, 6}. 
Then there exist positive constants so{m) and M — M{m) such that for s > so(m) we have: 

^1/2 f f,2s^m<^Q^z''{0) - dtZ''{0)f dx < M r I e^'^lg"" - g''\^ dxdt, (1.21) 
Jn Jo Jn 



where <p = (pi3,x is chosen so that Theorem \1.2\ holds and so(m) is the parameter in Theorem 

m 

Based on this result, we propose an algorithm to compute the potential Q from the mea- 
surement of the flux, based on some additional knowledge on the L°° (Sl)-norm of the unknown 
potential Q: 

3m > such that Q G L<,„(n). (1.22) 

The algorithm then is the following: 

Algorithm 1. • Initialization: q" — 0. 

• Iteration. Given , we set /i*" = dt {duw[q''] — duW[Q]) on To x (0,T), where w[q''] 
denotes the solution of 

( dtw- Aw + q''w = h, mnx{0,T), 

I w = ha, on9nx(0,r), (1.23) 

[ to(0) = Wo, dtw{0) = wi, in SI, 

corresponding to (|1.15p with the potential q^ . We then introduce the functional [/i* , 0] 
defined, for some s > that will he chosen independently of k, by 

Jj, gfc [/i*, 0](2) = I I e^^'^ldf z — Az + q'' zf dxdt 



Jn 

T 



+ \ [ [ e^'"^\d^z-t^''f dadt, (1.24) 
2 Jo Jrn 



on the trajectories z G L'^{0,T; IIo{Q)) such that d^z - Az + q''z G -L^(S7 x (0,r)), 
d^z G L^(ro X (0,r)) and z(-,0) = in Q. 

Let Z*" he the unique minimizer (see Proposition \4- -?[ ) of the functional qfc[/i'°,0], and 
then set 

~k+i^^k^dtZ%&l^ (1.25) 
wo 

where wo is the initial condition in (jl.lSp . 
Finally, set 

q'+'=T^iq'+'), wherer„(g) = | , f (1-26) 

^ ' ^ ' sign(g)m, if \q\ >m. 

One will see in Section |4] how Theorem 11.41 allows to prove the convergence of the above 
algorithm for s large enough: 

Theorem 1.5. Assuming the multiplier and time conditions (|1.2l) - (|1.3p and (|1.17|) . (|1.18|) and 

(|1.22ll . there exists a constant M > Q such that for all s > so(m) and fe G N, 



/ 

Jn 



In Jn 
In particular, when s is large enough, the above algorithm converges. 
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We emphasize here that this approach is constructive. In particular, at each step of this 
algorithm, we solve a quadratic strictly convex minimization problem, which can be easily 
done. Besides, the algorithm necessarily converges to Q due to Theorem 11.51 This is a great 
advantage compared to the classical methods for solving this inverse problem, that usually 
consists in minimizing 



w[q\ being the solution of p.l5|) (or p.23p ) corresponding to the potential q, which is not 
convex and may have several local minima. Of course, due to that fact, it is very difficult to 
propose a convergence result based on the minimization of the functional J in (|1.28|) . since 
classical minimization algorithms are not guaranteed to converge toward the global minimum 



However, when doing numerics, as underlined from the seminal work [15' to the most recent 
developments |13| , the discretization process usually creates spurious high-frequency waves that 
do not travel and strongly disturb control processes, even making the discrete controls blow 
up for some initial data to be controlled. Based on the recent work ^ that proves Carleman 
estimates for discrete waves uniformly with respect to the space discretization parameter, we 
will investigate the numerical methods to compute approximations of potentials in a future 
work. 

Let us finally conclude this section by giving some references considering inverse problems 
for hyperbolic equations using Carleman estimates. 

The use of Carleman estimates to prove uniqueness results in inverse problems was intro- 
duced in by Bukhgeim and Klibanov. The first proofs of the stability of inverse problems for 
hyperbolic equations rely on uniqueness results obtained by local Carleman estimates (see e.g. 
[171 122| ). One can read for instance [251 1261 [29l 1301 127] . where the method uses compactness- 
uniqueness arguments based on observability inequalities. 

Concerning other inverse problems for the wave equation with a single observation, the 
references [181 1191 120] consider the case of interior or Dirichlet boundary data observation and 
use global Carleman estimates, as in 2i. 

For more generic hyperbolic models, one could also mention [31 [S] [5T] giving stability of 
inverse problem from appropriate global Carleman estimates respectively for a network of 1-d 
strings, a discontinuous wave equation or the Lame system. Let us also mention the work [6] 
for logarithmic stability results when no geometric condition is fulfilled. 

1.4 Outline and notations 

The paper is organized as follows. Section [5] is devoted to the proof of several weighted esti- 
mates yielding the Carleman estimates of Theorems ll.ll and ll.2l In Section[3l we show how the 
result of Theorem 1 1 . 1 1 can be used to solve the exact controllability problem related to equa- 
tion (|1.10p . In particular, we give the proof of Theorem 11.31 on the dependence of the control 
with respect to the potential. Section [4] then focuses on the application of Theorem ll.2l to the 
inverse problem related to equation 1)1. 15[) of recovering the unknown potential. We describe 
in detail the algorithm proposed for solving this inverse problem and prove its convergence as 
stated by Theorem 11.51 

Due to the important number of notations, let us make precise some of them: 

• y represents controlled trajectories and u stands for the controls; 

• z represents free trajectories of the waves with a source term and satisfying homogeneous 
Dirichlet boundary conditions; 

• V represents trajectories of the waves satisfying v{±T) — dtv{±T) — 0; 

• w represents trajectories of the waves conjugated by the Carleman weight in Section [2l 
In the rest of the article, it represents solutions of 1)1. 23p : 




(1.28) 



of J. 
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• ip, tp are the Carleman weights defined by (|1.5P : 

• s,X,j3 are the parameters entering into the Carleman estimates; 

• p, q, Q denote potentials, but potentials denoted by p may depend on x and t whereas 
potentials q depend only on x. 

2 Weighted estimates 

In this section, we will prove several weighted estimates, including a weighted Poincare inequal- 
ity, in the goal of proving Theorems 11.11 and 11.21 

Let Q he a, smooth bounded domain of R", n > 1. In order to simplify the notations, we 
introduce the d'Alembertian operator: 

a = dt - A. 

In the following, we consider a function v € L^{-T, T; -ffo (fi)) such that Dv € L'^{Qx {-T, T)) 
and V = on d^l x {~T,T). We will first prove Carleman estimates similar to the ones of 
Theorems 1 1 . llfL2l for the operator □, corresponding to a potential equal to zero. 

2.1 A Carleman estimate in arbitrary time T 

Let us now give a global (meaning "up to the boundary") Carleman inequality, following 
Imanuvilov's method |17| . 

Theorem 2.1 (see [2]). Assume the Gamma- condition (|1.2p . Let ifj and he the weight 
functions defined by (II. 5p . 

Then for every (3 £ (0, 1), there exist Ao > 0, so > and a positive constant M such that 
for all s > So and A > Ao, 

sA r [ ipe^'^^ildtyf + \Vvf)dxdt + s''\^ f ip'^e^'^lvf dxdt 

+ [ [ \Pi{e'"^v)\'^ dxdt< M [ [ e^'"^\avf dxdt + Ms [ [ ipe^"'^ \dyv\^ dcrdt (2.1) 
J-TJn J-TJil J-tJto 

for every v G L"^ {~T,T; H^{D,)) satisfying Dv G L^(f^ x (-T,T)), d„v G L^(ro x (-r,r)) and 
v{±T) — dtv{±T) — m Q, and where Pi is defined by 

Piw = dtw ~Aw + s'^X'^ip'^w{\dtipf - 1 Vi/jj^). (2.2) 

Let us emphasize that such estimate is not new. Firstly, the proof can be read in the 
unpublished work [2, by the first author. Secondly, there are many Carleman estimates for 
hyperbolic equations. One can find (local) Carleman estimates for regular functions with 
compact support in [9l 1161 1141 1281 129) . One can also read similar versions of global Carleman 
estimates for hyperbolic equations in [171 1201 [3T] . 

For completeness, we give the proof of Theorem 12.11 below. 

Remark 2.2. This Carleman estimate is proved for any arbitrary time T > 0, but v has to 
satisfy v{±T) — dtv{iT) = 0. Therefore, the uniqueness result implied by Theorem \2.1\ is not 
surprising since the corresponding unique continuation result is: If v{±T) = dtv{±T) — 0, 
V G L'^{{~T,T)-H^{n)), Dv = and d^virg = 0, then v = 0. 

Proof. Using the weight ip defined by 1)1. 5[l . we set, for s > 0, 

w{x,t) = v{x,t)e'"^''''-*'> for all {x,t) e^lx {-T,T). 

Then, we introduce the conjugate operator P defined by 

Pw = e'^D{e~"^w). (2.3) 
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Some easy computations give 

Pw — dfW — 2s\ip{dtwdt'4' — Vw • V^) + w[\dt'(l^\^ ^ l^V'l^) ^ ^™ 

—sXip'w(dt^p — Alp) — sX^ipwQdt^l'l^ — |VV'|^) 
= Piw + P-zw + Rw, 

with Pi defined by and 

P2W = {a — l)sXipw{d1'4! — Aij)) ~ sX^ipw{\dt'>l>\^ — {Vipl^) 

-2sXip{dtwdtip - Vw ■ Vtp) (2.4) 
Rw = -asXifw{dttp - Alp), (2.5) 

where 

condition that will be explained later below in Step 2 of the proof. Since we have 

f f {\Piwf + \P2wf) dxdt + 2f f PiwP2wdxdt^ f f \Pw - Rwf dxdt, (2.7) 

the main part of the proof is then to bound from below the cross-term 

rT f 

P\wP2'W dxdt 



T Jil 

by positive and dominant terms, similar to the one of the left hand side of (|2.H) . and a negative 
boundary term, that will be moved to the right hand side of the estimate. For the sake of 
clarity, we will divide the proof in several steps. 

All the computations below are done for smooth functions v (equivalently, vj). Then, by a 
classical density argument, we can extend the results to any v satisfying Uv {-T,T)), 
d^v G L^{dn X (-r,T)) and v{±T) = dtv{±T) = in O. 

Step 1. Explicit calculations 

We set 

3 

{PlW,P2w)i2(nx{-T,T)) = ^'-^ 

where 7i,fe is the integral of the product of the ith-term in Piw and the fcth-term in P2'w. We 
mainly use integrations by parts and the properties of w such as w{±T) = 0, 9t«;(±T) = in 
n and w = on 9f7 X (-T, T). 

We shall also persistently use the fact that dt^ does not depend on x and that Vip does 
not depend on time, and thus 

dtV^ = dtVi^ = 0, 9tAV' = at(|ViA|') =0. 
Integrations by part in time give easily 



111 ~ I dtw{{a — l)sXipw{dt'4> — A^)) dxdt 
-T Jn 



= {l-a)sX r f 
J -T Jn 



ipldtwl^ [dtip — Aij)) dxdt 



~T Jn 

^^^"^ sA^ / [ ^IwfdtipidtTjj - AiP)dxdt 
2 J^TJn 



.(1^,,3 



X^ [ [ip\w\'^\dfipf{dtip-AiP)dxdt. 
J -T Jn 
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Similaxly, one has 



/i2 = / / diw{-s\''fiw{\dti)\^ - \^i)\^))dxdt 



/ / (fildtwi^ddtipf — \Vi!f)dxdt — sX^ / / (p\wf\dt^fdxdt 
J-tJu J-TJn 

-{2+hsX^ r f ^\wf\dtijfd^'iPdxdt+^ r f ^IwflV-ipfd^ipdxdt 

,4 f-T 



J !• 



' ' ip\w\'\dti^\\\dt^i>\' -\^M')dxdt 



and 



/i3 = / / dt 'w{—2sXip{dtwdttp — Vw • Vip)) dxdt 

-T Jit 

fT f fT 



-T Jn 

We compute in the same way 

I2I 

-T Jn 



sX [ I ip\dtw\^di'4>dxdt + sX" [ [ ip\dtwf\dtipf dxdt 
J —T Jn J —T Jn 

+sX f [(pldtwfAtpdxdt + sX^f f ifi\dtwf\Vtpf dxdt 
J-rJn J-rJn 

—2sX^ / / ifidtw dttp^w ■ dxdt. 
J-T Jn 



= f f -Aw{{a-l)sXipwidt^- A^))dxdt 
J-T Jn 

= -{l-a)sx[ [ ip\Vwf{d^tp-A'tP)dxdt 
J-T Jn 

_^(l_a)^_^2 /■ f ^\u,\^Ai>{dti> - ^Tp)dxdt 
^ J-rJn 

sX^ f f ip\wf\Vipf{dtip- AiP)dxdt 
J-T Jn 



and 



I22 = I I -Awi-sX^ifiwildttpf -\Vi;f ))dxdt 
-T Jn 

sX^ ( ( ip\Vw\'{\dt'il)\^ -\Vi)\^)dxdt 

J-T Jn 



_sX 
2 

sX 



T Jn 
2 rT 



^\wfA{\Vijf)dxdt 
T Jn 



f f ^\wfAiP{\dttpf-\V'ipf)dxdt 

J-T Jn 



T Jn 

4 rT 



+ ^ / / ^\wf\Vij\\\dti>f-\Vi;f)dxdt 



T J n 

T 



-sX^ [ [ (fi\wfVip-V{\Viljf)dxdt. 
J-T Jn 
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Using the fact that w|ar2x(-T,T) = 0, Vw = {dvw)v and |Vw|^ = on dQ. x {—T,T): 

I23 = f f -Aw{-2sX(p{dtwdtip - Vw • VV')) dxdt 

= s\ [ [ ip\Vw\\dtip - Ail))dxdt + 2sX^ f f ip\Vip ■ Vwf dxdt 
J-T Jn J-T J it 

-2sA^ / / (pdtwdtipVw ■ Vi^dxdt + sX^ I I ipl Vw|^(|9tVI^ - \^i'\^)dxdt 

-s\ f f (pld^wfVip-i^dadt + 'isX f f ip\Vwf dxdt 
J-T Jen J-T Jet 

since D^ip = 2Id. 
One easily writes 

hi = ( [ s'^X'^cp'^wildtipf -\Vipf){{a-l)sXcpw{dtip - AiP))dxdt 

J-T Jn 



-T Jn 

= {a-l)s^X^ f f (p^\w\\dtip - Aij){\dtijf -\Vijf)dxdt 
J-T Jn 



and 



I32 = [ [ s^X''^^w{\dti;f -\Vipf)i-sX^cpw{\dtipf -\Vipf))dxdt 

J-T Jn 

= -s^X" r I </\w\'{\dtil)\^ -{Vijff dxdt. 
J-T J a 

Finally, some integrations by part enable to obtain 



133 = I I s^X^-/w{\dt^\^ - |VV'|^)(-2sA<^(atw9tV' - Vw • ViA)) dxdt 
J-T Jn 

s^X^ r f ip^\wf{dfi;-AiP){\dtipf-\Vipf)dxdt 
J-T Jn 

+2s^X^ [ f ip^\w\\dtip\dtipf + 2\Vip\^)dxdt 
J-T Jn 

+3s^A* r f ip^\wf{\dtipf - \Vi;ffdxdt. 
J-T Jn 
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Gathering all the terms that have been computed, we get 
P1WP2W dxdt 



J -T Jn 



2s\ 



+ 2s\ 



+ 4sA 



' ip\dtw\^dttp dxdt — as\ / / ip\dtw\^ (d'ftp — A^) dxdt 
n J-TJn 

f f ip{\dtwf\dt^\^ -2dtwdt'ipVw-ViP + \V(P-Vwf) dxdt 

f f ip\\7wf dxdt + asX f f (p\\7w\'^ {dtip - Atp) dxdt 
J-TJn J-TJn 

fT 



- sA 



-T J an 

+ 2s^A^ r f ip^\wf{\dtipf - dxdt 
J-T Jn 

+ 2s^X^ [ [ ip^\wf{dtTp\dtipf + 2\ViPf)dxdt 
J-T Jn 

f f ip^'lwfid^i,- A^){\dtijf -\Vi:f)dxdt+ Xi, 

J -T Jn 



(2.8) 



+ QS^A^ 



/ / Lp^\wfdxdt (2.9) 
J ~T Jn 



-T Jn 

where Xi gathers the non-dominating terms and satisfies 

\Xi\ < MsX^ [ [ ip\wfdxdt. 
J ~T Jn 

Note that, since 1/' > 1, we have, for A large enough, A** < e^^^ — Lp^ , and therefore, 

/•T 

\Xi \ < Ms 

for some M independent of s and A. Here and in the rest of the proof of Theorem 12.11 M > 
corresponds to a generic constant depending at least on Q and T but independent of s and A. 

Step 2. Dominating terms 

On the one hand, about the first order derivative terms, one can notice that 

<fi {\dtwf\dtipf - 2dtwdtil^\/w ■ + ■ \/w\^) dxdt 

if {dtwdtip ~ Vw ■ Vtpf dxdt > 0. (2.10) 

T Jn 

Since this term can vanish, we focus on the terms in sA in l^tuij^ and | Vtyj^ and we want them 
to be strictly positive: 

2dfil) - a{dfip - Atp) > and 4 + a(9tV - ^■>P) > 0. 

By explicit computations, these two terms are positive if and only if a satisfies (|2.6p . This 
justifies the choice of the parameter a. Note also that this can be satisfied only if /3 £ (0, 1). 
As a direct consequence, we can write 



2sA 



// ipldtwl^dfip dxdt — asX / / ip\dtw\^ (df ip — A^p) dxdt 
-rJn J-rJn 

+ 4sX f f (p\Vwf dxdt + asX f f (f\Vwf {dttp - Aip) dxdt 
J-rJn J-rJn 

f ip\dtwf dxdt + MsX f f ip\Vwfdxdt. (2.11) 
in J-rJn 



> MsX 
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On the other hand, considering the 0-th order terms, we can observe that: 

T 

T Jn 



I 3x3 

+ as 



X' I I f>^\w\\d^t^ - A^)(|9tV!' - IVV!') dxdt 
J -T Jn 

2s^\^ [ [ ip^\w'^\{dfi>\dtipf + 2\Vipf)dxdt 
J -T Jn 

= s^A^ / / ip'^\w\^ F\{(t)) dxdt, 



where 

Fx{<l)) = 2A(19tVl' - IVVI')' + 2{dU\dti^? + 2|V^/>|') + a(9tV - A7/.)(|9tVl' - IV^/'l") 

= 2\{\dti>\^ - \'^i>?f + (29tV + a{dti> - Ai/>))(|9tiA|^ - |V7/>|') + 2(9tV + 2)\V^\' 
= 2AX^ + {2dt^ + a(9tV - A'0))X + 16(1 - li)\x - xof . 

with X = IVV'p. 

Since xq ^ Q, and /3 G (0,1), we have 16(1 — I3)\x — xa'\^ > c* > 0. Therefore, we are 
considering a polynomial A{X) > 2\X^ — 2 (q(/3 + n) + 2/?) X + c, and taking A > large 
enough, the minimum of A will be strictly positive. Consequently, 

s^A^ / / ip^\w^\Fx((l))dxdt> Ms^X-^ [ f i/\w\^ dxdt. (2.12) 



Thus, plugging (f2T0)) . ((2TTT1) and (t2T2)) in (p^Sjl . and since Vip = 2{x - Xo), we obtain 

/ P1WP2W dxdt + 2s\ / / ip \d,^w\'^ (x — Xo) ■ i'{x) dadt + \Xi\ 
'T Jn J-T Jan 

>MsX [ I if{\dtwf + \Vwf)dxdt + Ms^\^f I ip^\w\^dxdt. 
J-rJn J-rJn 



Since we also have 

f\Pw-Rwfdxdt < 2f f\Pwfdxdt + 2f [ \Rwf dxdt 
T in J-rJn J-rJn 

< M I I \Pwf dxdt + Ms^\^ [ [ if^\wfdxdt, 
J-rJn J-rJn 

using (|2.7p and (|2.9|l . we get 

s\ [ [ {\dtw\^ + \Vwf) ipdxdt + s^X^ [ f \w\^i/dxdt 

+ [ [ {\Piw\^ + \P2wf) dxdt 
J -T Jn 



<M f f \Pwf dxdt + MsX f I ip\d„wf (x - Xo) ■ u{x)dadt 
J-T Jn J -T Jto 

+ Ms f f ip^\wf dxdt + Ms^X^ f f ip'^\wf dxdt. 
J-rJn J-rJn 
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We take now so large enough so that the terms of the last line (coming from Xi and |-Rty|^) 
are absorbed by the dominant term in s''A^|w|'^<^'^ as soon as s > so- Using also the condition 
(II. 2|) on To, we finally obtain for some positive constant M, 

sX / if{\dt'w\^ + \\7'w\^)dxdt + s^X'^ / / ip'^l'wf dxdt 
J-TJsn J-TJsn 

+ [ [\Piwfdxdt+ [ f \P2wf dxdt (2.13) 



J -T Jn J -T Jn 

<Mj [ \Pwf dxdt + MsX f I ip\d^wf dadt 



T Jn 



for all s > So and A > Ao. 

Step 3. Back to the variable v 

Since w — we^*", we have 

e^^^ldtvf < 2\dtw\^ + 2s'\dtipf\wf < 2\dtw\^ + Cs' AVVi^ in Q x (-T, T), 
e^'"^\Vv\^ < 2\Vw\' + 2s^\V^\''\w\'' < 2\Vwf + Cs^X^ip'^\w\\ in x (-T, T), 
e^"^ \d^v\^ = \d^w\^ , on c»n X (-T,r). 

Using (|2.3p that gives by construction Pw — e^'^Ov, we can go back to the variable v in (|2.13|l 
and obtain that there exists some positive constant AI such that for all s > sq and A > Ao, 

r-T 



sX [ [ ipe^'"^{\dtv\^ + \Vv\^)dxdt + s''X'' [ f ip''e^'"^\v\^dxdt 
r I \Pi{e'^vfdxdt+ r f \P2{e'"^v)\'' 

I ipe'^'''' dadt. 



+ / / \P^{e"^ v)^ dxdt ^ \ \ \P2{e'^v]\'' dxdt 

rT r rT 

<M e''"^\Dv\^ dxdt + MsX 

J-TJn j-tJfq 

This concludes the proof of Theorem 12.11 □ 

In the sequel, we will fix A = Ao and use the fact that ip then is bounded from below by 1 
and from above by some constants depending on A. Since A is fixed, we can put it into the 
constants and obtain the following result: 

Corollary 2.3. Assume the Gamma-condition (|1.2|) . 

Then for every /3 G (0, 1), there exist A > 0, so > and a positive constant M such that for 
all s > So, 

s r f e'^'^ildivf + \Vvf)dxdt + s^ r f e''''^\vfdxdt+ I \Px{e'^v)\'' dxdt 
J-TJn J-TJn J-tJu 

<m[ I e^''^\nvf dxdt + Ms [ [ e^""^ \d^vf dadt (2.14) 
J-tJq J-tJto 

for every v G L'^ {{-T,T); Hi{Q.)) satisfying Dv G L^{n x {~T,T)), d^v G L^{dQ. x (-r,r)) 
and v{±T) — dtv(±T) — in Q,, where Pi is defined in (|2.2I) . 

In the following, A'l denotes various constants that do not depend on the parameter s. 

2.2 Weighted Poincare inequality 

We prove here a weighted version of the Poincare inequality that will be used thereafter. 
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Lemma 2.4. Let G C'^{Vi) and assume that the weight (p defined on Q, satisfies 

inf |V(^1 > 5 > 0. (2.15) 

Then there exist so > and M > such that, for all s > so and for all z £ _ff(5(f2), 

e^''^\zfdx<M f e'^'^lVzfdx. (2.16) 
! in 

Proof. We begin with the following computation, that uses the identity V • (e^'"^V(^) = 
e2''VA<^ + 2se2"^|V(^p: 

e''''^\zf\Vv>f d-x ^ ^J^ \z\^ (V ■ (e^'^^V^) - e'^^A^) dx 



= -s / e'"^zVz-Vifdx- - / e-'^'^lzl^Aifidx. (2.17) 
Since Lp £ C^{Q) and since we suppose (|2.15[l . the last term in (|2.17p can be bounded as follows: 

_£ f e'^"^\z\^Aifdx < Ms [ e'"^\z\^\V<^\^ dx. 



Then, for s sufficiently large, this term is absorbed in the left hand side of (|2.17p . This implies 
the following estimate: 



/ e^'^\zf\Viffdx<Ms [ le'^^zVi^l le'^'^Vzl dx 
in in ' Ml 



This yields (|2.16p and concludes the proof of Lemma [2.41 □ 

Let us emphasize that the weight function </p defined by (|1.5p is such that for all t G (— T, T), 
(^(•, t') satisfies assumption (|2.15p of Lemma [2. 41 since xo and Vi^ = 1\{x — xo)f. 

2.3 A Carleman estimate in time T large enough 



When the time T is large enough in the sense of p.3p . we claim that the conditions at times 
±T can be removed of the assumptions of Theorem 12.11 Roughly speaking, this will follow 
from an energy argument coupled to the Carleman estimate 1)2. ip . The result is given in the 
following theorem: 

Theorem 2.5. Assume the multiplier condition 1)1. 21) and the time condition (|1.3|) . Define the 
weight functions ip as in p.Sp with /3 G (0, 1) betng such that (|1.6|) holds. 

Then there exist sq > and a positive constant M such that for all s > so / 



e^'^ {\dtzf + \Vzf) dxdt + s^ r f e^'^lzfdxdt 

<M f f e'^'^\azf dxdt + Ms f f e'^'^ \d^z\'^ dadt, 
J-T in i-T irn 



(2.18) 



for all z G L"^ {{-T,T); Hl,{n)) satisfying Dz G L'^{n x {-T,T)) and d^z G L'^{dD. x {-T,T)). 



Remark 2.6. Let us emphasize that Theorem \2.5l contrary to Corollaru \2.3[ does not require 
z{±T) = dtz{±T) ^0 mn. 



15 



Proof. From condition (|1.6p that states sup^gQ \x — xo\ < /3T, we can choose rj £ (0, T) and 
e G (0, 1) such that 

(l-e)(r-77)/3> supjs-a;o|. (2.19) 
Then, explicit computations on tp{x,t) = \x — xop — jSt^ + Co show that we have 

C {l-e)\dtip{t)\>sup\V^{x,t)\, 
1 snpp(x,t) < Co < mi p(x,0). 

Hence, we introduce the cut-ofl function x G C*^ (R) such that < x 5- 1 a-nd 

" \ 0, if t<-T or t>T, ^^■^^> 

and we set v = x^, in x {—T,T). Therefore, v satisfies the required hypothesis v{±T) = 
dtv{±T) = in n and we can apply the Carleman estimate H2.14[l of Corollary 12.31 to v: 



s r I e^''^ {\dtv\' + \^v\'') dxdt + s^ r [ e^^'^lvf 
J-TJn J-TJn i-o ooN 

<M e^^'^\Uv\^ dxdt + Ms / e^'"^\d^v\^ dadt. 

J -T Jn J -T Jfo 

One can calculate that 

av = x^z - "ixdtz - x'z, in n X (-T, T). 

Besides, the functions x' and x" have compact support in (— T, — T + ry) U (T — 77, T). Thus, 
from (|2.22p . we deduce the following estimate on z: 

/ e'^'^ {\dtzf + \Vz\'^) dxdt + / / 
-T+T] Jn J -T+i-f Jn 

T 



< M 



[ [ e'^^lDzf dxdt + Ms [ [ e'^'^ld^zf dxdt 



+ M f [ e^'''' {\dtzf + \zf) dadt + M [ [ e'^'^ {\dtz\^ + \zf) dxdt. (2.23) 

JT-r] Jn J -T Jn 

We will show that the last two terms of (I2.23[) can be absorbed in the left hand side if the 
parameter s and the time T are chosen sufficiently large. In order to do that, we introduce the 
following weighted energy: 

Es{t) = i ^e^'^'*' {\dtz{t)\' + lVz(t)n dx. Vt e (-T,T), 

We first calculate: 

dE, 
dt 



^sl e^'^dfpildtzl^ + \Vzf)dx+ / e^''' {dtzdtZ + Vz-V{dtz))dx. 

Jn 

Thus, after an integration by parts, 

/ e^"^dtip{\dtz\^ + \Vz\^) dx + 2s [ e^"'^ dtzV ip ■ V z dx = [ e^'"^dtzDz dx. (2.24) 
dt Jn Jn Jn 

Step 1: Term of ^^23} on (T - ri,T) 

(a) Thanks to the formula 2ab < + b^, we can bound by below the left hand side of 
equality (|2.24|) leading to: 

^-s f e'^^^idt^+lVifDlldtzf + \Vzf)dx< [ e^'"^dtzUzdx. 
Jn Jn 
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According to (p:20)) and for te {T - r],T) 



int {-{dtip + \V If \)} > hif{-edtip} > 2e/?(r - 7])e^'^ > 2e/3(T - 77) := c. > 0. 



Thus, 



m€{T-r,,T), ^ + sc, f e^'"^ {\dtzf + \Vzf)dx< f e^'^^dtzDzdx. (2.25) 
™ Jsi Jn 

Now, using the formula 2ab < ea? H with e — sct, we can bound the right hand side of 

e 

(|2.25|l as follows: 



e'^"'^dtzOz dx 



2sc, 



The first term of the right hand side is absorbed by the left hand side of (I2.25|) . We obtain 

dE„ sc, .0 ,9x , , 1 

In 



or, equivalently, 

dEs „ 1 

— + sc^Es < 

dt Zsct 

Using the Gronwall lemma, we can write, for all t £ {T ~ rj,T), 



^sc.^r-t, , e^'"P(->\az{r)\UxdT 
'T-77 Jn 



Es{t) < Es{T-n)e'*^'^-^-'^ + f 

< i;,(r-r?)e-""*(*-(^-''» + -!- /' f e^'^^^\nz{T)\'dxdr. (2.26) 



Integrating this relation for t between T — rj and T, we obtain: 



I Es{t)dt < Es{T~rj) f e-=^*(*-(^'''''df + -^ / / 

JT-ri JT-7-1 2SC, J_yjQ 

< ^Es{T-r,) + — r f e^'^ldzf dxdt. (2.27) 
s s J-TJn 



(b) Now we want to estimate Eg (T — rj) by Es (r) for r G {—T + i],T — rj) . We use equality 
(|2.24[) that we integrate between r and T — rj: 



(T - ri)-Es(T) = s j " y e^^'^dtip {\dtz\' + IVz^ dxdt 
-2s / / e^'"^dtzVip-Vzdxdt + 



E, 



Jt Jn 

We bound the right hand side using Cauchy-Schwarz: 

r-T-v 



E,{T - r;) - E,{t) < Ms f " E,{t)dt + — f f e^'^\ 

J -T+ri J-TJn 

Integrating for r between —T + r; and T ~ rj, since s is large, 

Es{t)dt+— I e^'"^\Uz\^ dxdt. (2.28) 
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(c) Finally, thanks to (|2.27|) and (|2.28p . we deduce 

/ e^"^ {\dtz\' + \Vzf) dxdt 

T-ri Jn 

<M r f e^"^ {\dtzf + \Vz\'')dxdt+— r I e^'"^\Uz\' dxdt. (2.29) 
Moreover, combining it with the weighted Poincare estimate of Lemma [2.41 we get 

sf [ e^"^ {\dtz\^ + \Vzf + s'^\z\'^)dxdt< Ms [ Es{t)dt 

< Ms [ Es[t)dt + M f f e'^'"^\Uz\^ dxdt. (2.30) 

Step 2: Term of ((2:23|l on (-T, -T + 77) 

We want to obtain the same results as previously but on the interval (— T, —T-\-ri). In order 
to do that, one can introduce z{x, t) = z{x, —t) and apply the above estimates to z (we make 
the change of variable t — >■ —t). Thus, equations (|2.26|) - (|2.28[) coincide with the following ones: 



sc,{-T+ri-t) 



+ - I I e''^'^^^^\nz{T)fdxdT, (2.31) 



-T+17 



\ft£{-T,-T + ri), Es{t) < Es{~T + r])e 



/ Es(t)dt < —Es{-T + r,) + — e'^'^iazf dxdt. (2.32) 

J -T ^ ^ J ~T Jn 

fT-v »,r rT r 

Esi-T + T]) < Ms Es{t)dt+— e'^"^\azf dxdt. (2.33) 

J -T + r] ^ J -T Jn 

Combining (f232| with (fT33|l . we deduce 

r-T + r, r 

/ / e'"^ {\dtzf + \Vzf) dxdt 

J-T Jn 

<M r [ e^"^ {\dtzf + \Vz\'')dxdt+— r f e'^'^iazf dxdt. (2.34) 
J~T+nJn * J^rJn 



''^ " / e^'-^ {\dtz\^ + \Vz\^)dxdt+— r 
Besides, similarly to (|2.30|l . using Lemma [2.41 we have 

-T + ri 

-T Jn 



s 



r-T+v r 

/ / e^"^ {\dtzf + \\/z\^ + s^\z\^)dxdt 

J -T Jn 



<Ms f " Es{t)dt + M f f e^'^'p 
J -T+T] J -T Jn 



z\^ dxdt. (2.35) 



Step 3: Conclusion 

Using the power of s in the left hand side of 1)2. 23|) and estimates H2.29|l and (|2.34|l . taking 
a large enough, we obtain 



^ " / e'""'^ {\dtz\^ + \Vz\'' + s^l"") dxdt 

-T+r] Jn 

<M I f e^''^\Uz\' dxdt + Ms f f e'^'"^\d^z\'^ dadt. 
J-TJn J-tJto 

Using then estimates (|2.30l) and (|2.35p . we immediately deduce (|2.18p . □ 
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Remark 2.7. Note that, following carefully the above proof, the Carleman estimate (|2.18ll can 
actually be slightly improved into 

s r [ e""'^ {\dtzf + \Vzf) dxtit + s"" r [ e^'^lzl^ dxdt 

(2 36) 

<M ( [ e'^"^\az\^dxdt + Ms [ [ e'^"^x{-tf\d.4' dadt, 
where x *s the cut-off function defined in (|2.2ip . 

2.4 A Carleman estimate with pointwise term in time — T 

The proof of Theorem 12.51 easily gives furthermore an additional weighted estimate of the 
solution at time — T: 

Corollary 2.8. Under the conditions of Theorem \2. 51 we also have 

s[ e^"'^^-'^^ {\dtz{-T)\^ + \Vz{-T)f)dx + s^ f e^'"^^-'^^\z{-T)f dx 
Jn Jn 

<m[ [ e'^'"^\azf dxdt + Ms f [ e'^"''' \d^z\'^ dadt, (2.37) 

for all z G L'^{{-T,T);Hi{Q)) satisfying Eiz e L^(f^ x {-T,T)) and d^z G L^idQ. x {-T,T)). 

Proof. Using inequalities (|2.31|l at t — —T and (|2.33|) . we obtain an estimate on Es{~T). We 
then deduce (|2.37[l from the weighted Poincare inequality of Lemma [2.41 and Theorem 12.51 □ 

Let us now conclude with the proof of Theorem 11.11 that will be our main tool for the study 
in Section [3] on the design of a constructive process for building controls that depend weakly 
on the potentials. 

Proof of Theorem \l.l\ The Carleman estimate of Theorem 11.11 for the operator □ + p with 
p £ L^^{^1 X (— T, T)) is a direct consequence of Theorem 12.51 and CoroUarv 12.81 noticing that 

in n x{-T, T), 

\az\^ < 2\az + pz\'^ + 2||p||ioc(ox{-T,T))|z|^ < 2\Dz+pz\^ + 2rn\zf. 
Then choosing so large enough, one can absorb the term 



e '"'^\z\ dxdt 



by the left hand side of the sum of (|2.18p and (|2.37p . thus obtaining (|1.7p with slightly different 
constants. □ 



Remark 2.9. Using Remark \2. 71 one can easily adapt the above proof to estimate the left 
hand-side of (|1.7[l by 

M f f e^'"^\Uz+pz\^ dxdt + Ms f f e^'"^x{'tf\d^A^ dodt, (2.38) 
J-tJsi J-tJto 

where x *s the cut-off function defined in (|2.2ip . 
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2.5 A Carleman estimate with pointwise term in time 

We are now interested in deriving a Carleman-type estimate with dtz{-, 0) in the left-hand side 
under the condition that z{x, 0) = for all x G Q. We aim at proving Theorem [L2] at the end. 

Theorem 2.10. Assume the multiplier condition (|1.2p and the time condition (|1.3p . Define 
the weight functions </p as in (II. 5|) with (3 £ (0, 1) being such that 1)1. 6|) holds. 

Then there exist so > 0, A > and a positive constant M such that for all s > sq: 

^1/2 f e2s^{0)|9j^(o)|2da; + s f e^'^ {\dtz\^ + \V zf) dxdt + s"" C f e^'^^izf dxdt 
Jn J-TJn i-T in 



T 



<M e'""^\az\" dxdt + Ms / e""^ dcrdt (2.39) 

i-T in J-TJTo 

for all z e L2((-r,T);//(l(n)) satisfying Uz G L^(f7 X {-T,T)), d^z G L^idQ. x {-T,T)) and 
z{x, 0) = for all x £ Q. 

Remark 2.11. Let us emphasize the assumption z{-,0) — in Q,. Without this condition, we 
do not know how to estimate the weighted energy Es at time t = 0. Whether this is a purely 
technical artifact or not is an open problem. 

Proof We consider a fmiction z G L'^ {{-T,T); Hl^{Q)) such that z{x,0) = for all x & D,. We 
use the notations previously introduced by (|2.2p and (|2.2H) and set 

w = e"'^x^ S'lid Piw = dfW — Aw + s^ X^(fi^w{\dttp\'^ — IVV'I'^). 

Since we are interested only in the dependence of s, as before, all the powers of A and the 
functions (p can be omitted and enter into the constants. 

Under the condition z{-, 0) = in Q,, we get w{x, 0) = for all x & Q. This allows us to do 
the following computations 

Or f'^ f 

P\W dtw dxdt ~ I I {d^w — Aw + s^ Lp^w{\dt'4'\^ ^ dtW dxdt 



T Jn 



\ £ \^^wm^ dx-'^ £^ j^^ \w\'dt {vWdU? - IVv^n) dxdt 



-T Jn 
rO 

> I / \dtw{0)f dx - Cxs^ 



1 



In 

implying in particular, by Cauchy-Schwarz, that 

s^^^ [ \dtw{0)\'^ dx < f f \Piwfdxdt + s f f \dtw\^ dxdt + Cxs^''^ I f \wf dxdt. 
Jn i-T in i-T in i-T in 

Moreover, v — lue^^*' = xz satisfies the assumption of Theorem 12. II Therefore we can use 
estimate (|2.13p on w and, bounding each from above and from below, we get: 

s^^^ f \dtw{0)f dx + s f f {\dtw\'^ + \Vwf) dxdt + s^ 
Jn i-T in 



T Jn 

+ [ f\Piwfdxdt+f [ \P2wf dxdt (2.40) 
i-T in i-T in 

<M f f \Pw\^ dxdt + Ms ( f \d^wf dadt. 

i-T in J-TJrn 



Now, arguing as in Step 3 of the proof of Theorem 12.11 we obtain 

gi/2 f e2-»'(0)|a^^(o)i'' ds + s r [ e''"^{\dt{xz)f + \V(xz)f)dxdt + s^ f [ e'""^x^\z\' dxdt 
Jn i-T in i-T in 

<M [ [e'^''-^\D{xz)fdxdt + Ms[ [ e^"^ x^ \d.zf dadt. 
i-T in J-tJto 

We then use energy estimates as in Theorem 12.51 to deduce (|2.39p . □ 



20 



Let us now conclude with the proof of Theorem II. 2 1 main tool to study the inverse problem 
in Sectional 

Proof of Theorem W^ Since p G i<„(r2 x (— T, T)), thanks to Theorem 12.101 and arguing as 
in the proof of Theorem ll.il the potential in (|1.8p can be absorbed by taking s large enough. 

When the potential q in the operator does not depend on time, one can extend the function 
z by z{-,t) = z{-,~t) for t £ (— T, 0) and apply (11.81) to this extended function z. Of course, 
since each term is odd or even, the integrals on (— T, T) simply are twice the integrals on (0, T), 
which concludes the proof of (|1.9p . □ 

Remark 2.12. Using Remark \2. 71 similarly as in Remark \2.9[ one can estimate the left hand 
side of 113} by (fOSl) . 



3 Application to a controllability problem 

In this section, our goal is to present what are the consequences of the Carleman estimate of 
Theorem 11.11 with respect to the control properties of equation (|1.10|l . 

In all this section, we shall assume that conditions (|1.2|| - (|1.3p on To, T hold. Then there 
exists P £ (0, 1) such that (|1.6|l holds. We fix /3 this way and take A large enough so that 
Theorem 11.11 applies. 



3.1 Setting 

Let us recall that the exact controllability problem under consideration is the one described in 
the introduction by l|1.10|l - (|l.ll[l . In order to solve that problem, following the duality tech- 
nique introduced in T4 for parabolic equations (and that can be seen as an extension of the 
usual Hilbert Uniqueness Method 124j), the idea is to minimize the functional Ks^p defined by 

Before going further, let us take some time to describe the space on which Ks,p is defined. 
In the introduction, for p G L°°{yi x (— T, T)), we defined Ka.p on 

T[p] = [z£ L\-T,T;H^{n)), with {d't -A+p)z e L'{Q x (-T,T)) 

and d^z G L'^iTo x (-r,r))|. (3.1) 

Note that this is a space of trajectories of the wave operator with potential {df — A + p) and 
therefore it a priori depends on p. In order to study the functional Ks^p, natural semi-norms 
on T[p] are the following ones: 



l^llobs,s,p = — / / e'""^\dt z — Az + pz]'' dxdt + I I e'""^\d,jz\'' dadt. (3.2) 
" -T Jn J-T Jrn 



s 



Let us explain below that these quantities define norms for any parameter s > and potential 
p G L°°{Q X {—T,T)) and give some of their basic properties: 

• For aU s > and p G L°°{^1 x (— T, T)), and for aU z G T[p], the quantity ||2||obs,s,p is 
equivalent to the quantity 

\\z\\lbs,p= [ f\dtZ-Az+pzfdxdt+f f \duz\^dGdt, (3.3) 

in the sense that there exists Cs such that for all z G T[p], 

'T^\\^\\ohs,p ll^llobs.s.p ^ C'sll^llobs.p- (3.4) 

This is a consequence of the fact that the weight function e^"'^ is bounded from below 
and from above by positive constants depending on s. 
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• For p G L°°{Q. x (-T,r)) with s > so(||p||L=°(nx{-T,T))) given by Theorem [LH the 
quantity defined in (|3.2p is indeed a norm: the Carleman estimate p.7p shows that if 
ll^llobs, s,p = 0, then z = and that ||z||obs.s.p measures the L'^(— T, T; //q (f2))-norm of z. 
This proves that || ■ ||obs,s,p is a norm on T[p] for all s > so(|lj3|lL={nx(-T.T)))- According 
to the first item, this is actually true for all s > 0. 

• T[p] = T[0] for allp e L°°(f2x(-r,r)). This is due to the fact that pz G L'^ {^Ix {-T,T)) 
whenp G L°°(D. x {-T,T)) and z G L^(-T, T; //(J(n)). 

For convenience, we shall now denote T[0] simply by T, where 

r = |z G L^{-T,T;HI{Q.)), with {di - A)z G L\n x {-T,T)) 

and G L^(ro x (-r,r))|. (3.5) 

• For all m > 0, there exist constants M(m) > and S()(m) > independent of s and p 
such that for aU s > so(m) and p'',p'' e -L<„(ri x (-T,T)), 

^INIIobs,s,p'> ^ ll^llobs.s.po < A/||2:||j,j,g (3.6) 

This result follows immediately from Theorem 11.11 and its proof. Note that in (|3.6p . these 
equivalences of norms are proven uniformly with respect to s > so {m) for potentials lying 
in L^„(S7 x {—T,T)). This is an important remark. 



In the following, we will denote by (T, || • ||obs,a,p) the space T endowed with the norm 
II ■ llobs.s.p. But we drive the attention of the reader to the fact that these norms || • ||obs.s.p are 
not uniformly equivalent with respect to s > 0. 

We now introduce the functional space Hq{Q,) x L^{Q,) on which we consider the family of 
norms, 

11(^0, ^i)||-T,. = / e^'''^-^^\Vzo\^ + \z^\^)dx, fors>0. 
in 

Again, several remarks can be done: 

• Using the weighted Poincare inequality of Lemma 12.41 || ■ \\~t,s is equivalent, uniformly 
in s > 0, to 

( f e''^^-^^{\Vzo\' + \zr\')dx + s' f e''^^-^^\zo\' dx 
\Jn Jn 

• According to Theorem [TTl for aU 2 G T, (z{-T),dtz{-T)) belongs to H^{Q.) x L^(n) 
and for all m > 0, there exists a constant C independent of p and s such that for all 
p G X (-T,r)) and s > so(m), 

\\{z{-T),dtz{-T))\\-T.s<C\\z\\,t„,s,p, zeT. (3.7) 

Note that, anyway, for all s > and p G L°°{^ x {—T,T)), one easily checks that, bounding 
the functions depending on s if needed, there exists C{s,p) such that 

\\{z{-T),dtz{-T))\\-T,s <C{s,p)\\z\\obs,s,p, zeT. (3.8) 

Finally, we also introduce the dual space L^(f2) x H^^{Q,) on which we consider the family 
of norms 

\\{yo,yi)f-T,s,,= |^e-^^^(-^)(lyo|' + |V(-A,)-V|')d^. (3.9) 
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3.2 Construction of a null-controlled trajectory 

Proposition 3.1. Assume the multiplier condition (|1.2p and the time condition (|1.6[l . 

Then, for all s > and p G L°°{Q x {—T,T)), the functional Ks,p defined in (|1.12l) . that 
we recall below for the convenience of the reader, 

Ks,p(z) — I I e'^^^ \dt z — Az + pz\'^ dxdt + ^ I I e^"'^\dvz\^ dadt 

+ (^(-r),at2:(-T)))(^2x/f-i)x(Hixi2), 

is continuous, strictly convex and coercive on (T, || ■ ||obs,s,p) for initial data {vq^ ,yi^) in 
(L^(n) X H~^{Q), ||(-, •)||_T,s,*) and therefore admits a unique minimizer Z[s,p] G T. 
Setting 

Y[s,p] = h'^'^idf -A + p)Z[s,p] and U[s,p] = e^'"^d^Z[s,p]lro 

as in H1.13p . solves p.lOp with control U[s,p] and satisfies the control requirement (|l.lip . 

Besides, for all m > 0, there exists a constant M > independent of s and p such that 
for all p G L<^(r2 x {-T,T)), s > so(m), for all data (y,7^,yr^) G L^{n) x H~'^{n), the 
minimizer Z[s,p] of K^^p satisfies: 

ll^[s,p]llob=.s.p = s / e-^''^\Y[s,p\\''dxdt+ r f e-^'^\U[s,p]fdxdt 

<M\\iyo^,yT^)f_r.,,,. (3.10) 

Proof We fix s > and p G x {-T,T)). For {yo'^,yi'^) G L'^{Q.) x H-\n), the 

functional Ks,p is defined and continuous on (T, || • ||obs,s,p) because of (|3.8p . Estimate (|3.8p 
also yields immediately the coercivity of Ks.p: 

Ks.p{z) > ^ii2:||obs,s,p - 11(^(7^, ^r^)ii-T,sii(y(r^,yr^)ii-T,s,* 

> ^PI|obs,s,p - C'(s,p)||2||obs,s,p||(y(7^,yr^)l|-T,s,*. 

Therefore /S's^p has a unique minimizer Z[s,p] on (T, || • ||obs,s,p) • 

Since Ka^p{0) = 0, we have Ks,p(Z) < 0, which, according to the above inequality, imphes 

\\Z[s,p]\\obs,s,p < 2C{s,p)\\{yo'^ ,y^'^)\\^T,s,*- 

Besides, when for some m > 0, we have s > so(m) and p G L^^{i^ x (— T, T)), using 
(|3.7p instead of H3.8p . the above constant C can be chosen independently of s > so(m) and 
p e L^^n X (-T,T)): 

\\Z[s,p]U,,,,p < 2C\\{yo^,y^^)\\-T,s,„ 

which is precisely p.lOp . 

Let us now check that F[s,p] defined by p.l3p is a controlled trajectory of p.lOp with 
control function (7[s,p] that satisfies the control requirement (jLlip . In order to simplify the 
notations, until the end of the proof, we fix s > and p G L°°{il x (— T, T)) and denote 
Z[s,p],Y[s,p],U[s,p] hyZ,Y,U. 

The Euler-Lagrange equation given by the minimization of Ks,p is as follows: for all z £ T, 

- f f e^"^{dtZ-Az + pz){dtZ-AZ + pZ)dxdt+ f f e^^^d^zd^Zdodt 
s J _T Jn J ~T Jtq 

+ {(y(7'^,yr^), (2(-r),9t2(-r)))(^2xH-i)x(HixL2) =o. (3.11) 
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Therefore, with Y and U as in (|1.13|l . we obtain, for all z £T. 

I {dtz - Az+ pz)Y dxdt + ( 
-TJn J~tJt, 



I I {dtZ - + pz)Y dxdt + f I duzU dadt 

+ {{yo'^,yi^), (2(-r),9t2(-r)))(L2xif-i)x(HixL2) = o. 



But this is precisely the dual formulation of equation H1.10|) and integrations by parts yield, 
for all z G T, 

/ f z{df - A + p)Y dxdt + [ f duz{Ulva-Y)d(jdt 
J -TJn J -T J an 

This implies that Y = Y[s,p\ solves (fLTO|) - (frTT|) with control function U = U[s,p\. □ 

3.3 Dependence of the controls with respect to the potentials 

We are now in position to prove Theorem 11.31 



Proof of Theorem \1.3\ In order to simplify the notations, we set — dt — A+p^ , — Z[s,p''], 
y = y[s,p'] and [/* = U[s,p^] for i £ {a.,b}. Using (|1.13|l . we are going to bound the following 
expression: 



T 

s I I e-^'"^\Y'' -¥"'['' dxdt + 
-T Jn 



I I e~'""^\U'' -U''fdxdt 

- r [ e^'"^\L''Z'' -L''Z''fdxdt+ r f e^"^\d,Z'' - d^Z''\' dxdt 
■5 J -T Jn J -T Jth 



TJn 

^ [ e""^ (iL^Z^f + \L''Z''\'' ~2L''Z''L''Z'') dxdt 
-TJn ^ ' 

T 



+ f I e^"^ ({d^Z'^l^ + \d^Z^\'' -2d^Z''d^Z^'\dxdt. (3.12) 

J — T J Fq 

We recall first the Euler-Lagrange equation (|3.1ip associated to the trajectory Z\ for i £ {a, fe}: 



e^'"^VzVZ' 
s J_T Jn J -T Jto 



dxdt + / / e^"'^ dvzduZ^ dxdt 
J -T Jtq 

+ {{yo^,yi'^),{z{~T),dtz{-T)))i^^2.^H-i)x(HlxL^) =0 (3-13) 



We apply (|3.13|) to z = Z^ and z = Z\ j ^ i and obtain the following equations for € 
{a,fef: 

e'''"^\VZ'fdxdt+ r I e^'^^ld^Z'l^dxdt 
TJn J-tJfq 

+ ((j/o~^,yr^),(^'(-r),a,Z'(-T)))(^.,^_i),(^i,^.)=0; (3.14) 



s 



f f e^"'^ Z' V Z' dxdt + I [ e'""^d^Z'd,,Z'dxdt 
J -TJn J-tJto 

+ {(2/o"^,yr^),(^'(-r),aiZ^(-r)))(^.^^_i)^(^i,^2) = 0. (3.15) 
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Summing (|3.14[) for i — a and i = h and subtracting from it (|3.15p for — {a,b) and 

= (6, a), we get: 



T Jn 



^.sv ( jL^Z"!" + \L''Z°Y - VZ^VZ'' - L" Z" L" Z" j dxdt 

+ r f e'"^ (jd^Z^f + \d^Z''f - 2d^Z''d,Z'') dxdt = 0. 
Thus, making use of UZ^ = U Z^ + {p^ -p')Z^, we have from ((3A2|) : 

f e-2=^|F"-y*|2d^dt+ f e-^^^^lU" ^U''f dxdt 
-TJn J-tJtq 



fT 

e^"'^ {L" Z'' Z" + L^'Z^L^'Z'' - 21" Z" L'' Z'') dxdt 

T Jn 

T 

' e'^'^ip" - p''){Z'' L" Z" - Z'^VZ") dxdt. 
T Jn 



Finally, we obtain 



e-^^^\Y--Y''\^dxdt+ I I e-"=^|(/"-C/''l^d:rdf< ^||p"-p°||ioc(nx(-T,T)) 

T Jn J -T J To s 

T r f-T . 



e''"^(\L''Z''\^ + \L''Z''\^)dxdt + s^ / e''"^(\Z"\^ + \Z''\'')dxdt] . (3.16) 
TJn j-tjtq ) 

Applying the Carleman estimate of Theorem II. II for i G {a, &}, we have: 

e^'''^\Z'f dxdt<M r [ e^'"^\VZ'\'' dxdt + Ms f I e^"^\d,Z'\'' dadt. 
TJn J -TJn j-tjtq 

Of course, this implies 

e^"^ 1 1 ' dxdt + S-' I I e"'^ I dxdt 
-T Jn 



<M f [ e^"^\L'Z'\' dxdt + Ms [ f e'"^\d,Z*f dadt 
J -TJn j-tjtq 



and since 



e^-^\L-Z'\' dxdt + s / / e^'''^\d^Z'\^ dadt 

TJn J -T JTa 

e-'^'"^\Y'\'' dxdt + s r [ e-^'^^lU'-f dadt, 
TJn j-tjtq 



it can be rewritten as 



/ f e^'^lVZ'l^dxdt + s^ f [ e^^'^lZ'fdxdt 
J -TJn J -TJn 



-TJn J-TJn 

f-T 



<Msl^sJ J e-^'^\Y'\Uxdt + J J e-^'^\U'rdxdtj. (3.17) 
Combining (|3.16p and (|3.17|l . we obtain the desired estimate p.l4p . □ 
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4 Application to an inverse problem 



Let us consider the inverse problem defined in (|1.15p . In this section, we shall propose an 
algorithm based on the Carleman estimate (|1.9p and a data assimilation approach. Let us recall 
that the unknown is the potential Q = Q{x), that we aim at recovering from the measurement 
of the normal derivative of the solution W^[(5] of 1)1. 15|) on To x (0,r). We also assume that 
Q £ L<'^(0) for some given constant m > 0. 

Let us also mention that we are working under the geometrical assumption (ll.2l) - (|L3|) . and 
the function ip we shall consider below always satisfies p.6p such that Theorem 11.21 holds when 
the parameter s is large enough. 

As said in the introduction, one can find in [2] the proof of the fact that the additional 
information 9yW^[Q] on To x (0, T) allows to identify Q uniquely within the class of potentials 
in L'^^{Q). Our approach will go further, providing an explicit algorithm to compute Q. 

Our goal is indeed to prove that Algorithm [1] presented in Section [LS] is convergent when 
s is large enough, as described in Theorem 11.51 

In the following, we shall first present the idea underlying this algorithm. We will then 
focus on the proof of Theorem 11.41 which is the main step within the proof of the convergence 
result of Theorem 11.51 

4.1 The general idea 

Algorithm [T] is based on the fact that if W is the solution of equation p.l5p and wlq^] solves 
(fOal) . then 

z'^ = dt [wlq"] - W[Q]) (4.1) 

solves 

^Iz*^ - Az*^ + g^z'-' = 5^ innx(0,T), 

z'' = 0, on dn X (0,T), (4.2) 

2*^(0) = 0, dtz'^iO) = zt inn, 

where 

g" = (0 - q'')dtW[Q], zi = {Q- q'')wo, (4.3) 

and by definition, 

^'^ = a^z*^ on To X (0,T). (4.4) 

Of course, both variables g'' and zf are unknown, but the variable g'^ brings lower order 
information than /j,''. This fact actually is the milestone of the proof of Theorem 1 1.5 1 We shall 
then try to approximate zf through the additional information (|4.4p and let the source term 
free, as it is in the functional [^*, 0] in p.24p . 

Note that this idea is behind the proofs of stability by compactness uniqueness arguments 
as in [251 [261 [M EZ] or by Carleman estimates given in [21 [TSl [T9ll20] . 

4.2 Study of the functional Js^q[fi,g] 

We first give a functional setting for the minimization of the functional Js,q [fi, g] given in (|1.20p . 
that we recall below for the convenience of the reader: 

Js,q\p, g]iz) = ^ / / e''^''^\dtz - Az + qz - g\'^ dxdt + I- / / e^^'^|9^z - dcrdt, 
Jo Jn ^ Jo Jvo 

defined on the trajectories z such that z G L^(0, T; (H)), 9|z - Az + qz £ L^(n x (0,T)), 
9^z G L^{Vo X (0,r)) and z(-,0) = in 

Of course, it is very close to the functional Ks^p and we shall therefore introduce the space 

r+ = |z G L''{Q,T;HI{Q.)), with a^^z - Az G L^(n X (o,r)), 

z(-,0) = 0inJ7 and a^z G L^(ro X (0,T))|, (4.5) 
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and the family of norms 



||2 



- [ [ e^''^\d^z- Az + qz\^dxdt+ [ [ e^'"^\d,z\^ dadt. 
* JQ Jn Jo Jtq 



Note that these are norms for all s > and q G L°°(S1) according to Theorem 11.21 

The properties of the space T'^ and the family of norms || ■ ||obs+,s,q ^ire of course completely 
similar to the ones of T in (|3.5[l endowed with the family of norms || ■ ||obs,s,p introduced in (|3.2p . 
Therefore, we refer the reader to Section |3] for remarks and comments on (7"^, || ■ ||obs+ s q)- 

The first result states the well-posedness of the minimization problem of Js,q[jJ.,g]. 

Proposition 4.1. Assume the multiplier condition (|1.2p and the time condition (|1.6|) . Assume 
also fi e L^(ro X (0,T)) and g £ L^(D. x (O.T)). 

Then, for all s > and q £ L°°{n), the functional Js,q[fJ,,g] defined in (ll.20p is continuous, 
strictly convex and coercive on {T'^,\\ ■ Hobs+.s.g)- The functional Js,q[fi,g] therefore admits a 
unique minimizer Z in . 

Besides, for all data {ii,g) G L^{dQ x (0, T)) x L^(n x (0, T)), the minimizer Z of Js,q [/i, g\ 
satisfies: 

< - r / e'-^\gfdxdt + A f e'^^M^ dadt. 



s 



Jr 



Proof. The continuity, strict convexity and coercivity of the functional Js,p[lJ.,g] is straightfor- 
ward and left to the reader. 

To get estimates on the minimizer Z, we use Js.q[fJ,, g]{Z) < Js^q[fi, g]{0): 



( e^'^^ldtZ - AZ + qZ - g\^ dxdt+ ( f e^'^^ldr^Z - fif dadt 
Jn Jo Jfq 



<- f f e^''^\gfdxdt+ f [ 
Jo Jn Jo Jrg 



T 

e '"^ dadt. 



Developing the square on the left hand side, we obtain 



e'^'"^\dtZ~AZ + qZfdxdt+f f e^'"''\d„Zf dadt 
Jn Jo J To 

<- f f e^'^^idtZ - AZ + qZ)gdxdt + 2 f I e'^'"^d,,Zn dadt. 



s 



Jn Jo Jr 



Using 2ab < + 26^, we thus obtain 



e'""^\gfdxdt + 2 / e^""^ \^\' dadt, 
Jo Jfo 



which concludes the proof of Proposition 14.11 □ 

Of course, our goal is not only to prove that the functional Ja,q[lJ.,g] has a minimum, but 
rather to study how the minimum of Js,q[n, g] depends on the source term g. Indeed, z'' in (|4.1[) 
is the minimum of the functional [fj,'' ,g''], whatever s > is, whereas in the algorithm, Z'' 
is the minimizer of the functional J^^^k [/i'' , 0], see (|1.24p . 

This is precisely the goal of Theorem 11.41 As in Section 13.31 we shall rely on the Euler 
Lagrange equations satisfied by the minimum of the functionals Js,q{fJ., g"'] and Js,q[fi, g*"]. 

Proof of Theorem \1.4\ Let us write the Euler Lagrange equations satisfied by Z-' , for j £ {a, b}: 

- r I e'^'^^idtZ' - AZ' + qZ' - g'){dtz - Az + qz) dxdt 
s Jo Jn 

+ [ I e'^^^id^Z-' - fi)d^zdadt = 0, (4.6) 
Jo Jfo 
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for all z £ T^. Applying (|4.6p for j = a and j = b to z = Z"^ ~ and subtracting the two 
identities, we obtain: 



/ e'^'"^\dtz- Az + qz\'^dxdt+ [ [ e'^'"^\d^zf dadt 
Jn Jo Jto 



— - / e^'"^{g'' — g''){dtz — Az + qz)dxdt. 

^ Jo Jn 

This implies that 



e''"'^\dt z ~ Az + qz\^ dxdt + s j / e^"'^\d^zf dadt 

Jfo 

< 



Uo l^^'^^\9'' - a'l' dxdt. (4.7) 



But the left hand side of 1)4. 71) precisely is the right hand side of the Carleman estimate p.9p . 
Hence, applying Theorem ll.2l to z, we immediately deduce (|1.2ip . □ 

4.3 Convergence of Algorithm 1 

Let us now focus on the proof of Theorem 11.51 



Proof of Theorem \1.5l We use Theorem 11.41 since, as we explained, we have to compare the 
minimum Z'' of J^^fc[/x'°,0] with z*" = dt{w[q'^] — W[Q]) solution of (14.21) . which corresponds 
to the minimum of p*]. Note that this requires q'^ £ L'^^{Q,), which is guaranteed at 

each step of the algorithm by (ll.26|) . We obtain 

^1/2 f e''''^(o)\g^z''{0)~dtz''{0)fdx<Mr[ e'""^\g'=\' dxdt. (4.8) 
in Jo Jn 

But, from (fOSj) and (fO)) . 

9tZ'=(-,0) = (g'+^~g")u;o, dtz\- ,0) = {Q ^ q'')wo and g" ^ {Q - q'')dtW[Q]. 
Therefore, since ip{-,t) < ^{-,0) for all t G (0,T), estimate (|4.8[) reads: 

Jn Jn 
Of course, using the strict positivity (|1.18|) of wo, this yields in particular that 



° ' (mfn Iwol) 



n 



Since Tm defined in p.26p is Lipschitz and Tm{Q) = Q (because Q £ L^^{Q)), we have 
|g''+^ - Qi = |T™(g'=+^) - T{Q)\ < \q''+^ - Q\, from which we immediately deduce ([OT)) and 
conclude to the convergence of Algorithm 1 for s large enough. □ 

5 Conclusion 

As a conclusion, let us formulate a few comments and highlight some remaining open problems. 

On the geometrical conditions. Our strategy requires the use of Carleman estimates, 
and in particular the conditions 1)1. 2p and p.3p . But these conditions are much stronger than 
the classical Geometric Control Condition (GCC) introduced in [T]. Whether or not similar 
results as the ones presented above apply when only the GCC holds is an open problem. In 
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particular, to our knowledge, the only stability result in inverse problem proved using micro 
local analysis is the recent work |27| . which requires the GCC and the convexity of the whole 
boundary. 

Smoothness of Controls. The control process proposed in Section [13] does not fit in the 
framework developed in 12 which proves that using the Hilbert Uniqueness Method (HUM) 
(slightly modified by the introduction of a smooth cut-off function in time) to compute the 
controls, if the data to be controlled is smooth, then the corresponding control and controlled 
trajectory are smooth. Therefore, new questions arise: 

• Does the control process in Section [L2] enjoy smoothness properties similar to the ones of 
the classical HUM control? Note that these regularity properties arise naturally when consider- 
ing the control properties of semi-linear wave equations - see |10| - or when deriving convergence 
rates for the discrete controls, as explained in |13| . 

• How does the usual HUM control process depend on the potentials of the wave equation? 

Numerics and inverse problems. Recently, in [4 , we have proved discrete Carleman es- 
timates for the space semi-discrete 1-d wave equation discretized using finite differences. There, 
following the results on the observability of discrete waves - see e.g. |13) -, a new term has been 
added to make the Carleman estimates uniform with respect to the discretization parameter. 
This term, somehow corresponding to some kind of Tychonoff regularization of the Carleman 
estimates, is needed due to spurious waves created by the discretization process. Based on 
these uniform Carleman estimates, we have been able to prove a convergence result for the 
approximation of a potential in the inverse problem given in Section [T31 provided a Tychonoff 
regularization term is added in the process. 

It would then be completely natural to try to adapt the algorithm developed here in the 
continuous case to the space semi-discrete schemes and in numerics. This is currently under 
investigation. 
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